Exact solution for two interacting electrons on artificial atoms and molecules in solids 
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We present a general scheme for finding the exact eigen- 
states of two electrons, with on-site repulsive potentials {Ui}, 
on I impurities in a macroscopic crystal. The model describes 
impurities in doped semiconductors and artificial "molecules" 
in quantum dots. For quantum dots, the energy cost for 
adding two electrons is bounded by the single-electron spec- 
trum, and does not diverge when Ui — > oo, implying limita- 
tions on the validity of the Coulomb blockade picture. Ana- 
lytic applications on a one-dimensional chain yield quantum 
derealization and magnetic transitions. 

31.25.Nj, 71.70.d, 73.20.Hb, 73.23.Ps 



There has been much recent interest in the effects of 
interactions on the localization of electrons in disordered 
systems An important question, triggered by recent 
experiments in two dimensions Q , concerns the possibil- 
ity that the interactions help to delocalize the electrons, 
yielding a metal-insulator transition. For two interact- 
ing electrons in a random potential, this possibility was 
supported by Shepelyansky's numerical work by a 
Thouless type block-scaling argument El and by further 
numerical scaling results (|||. The latter took advan- 
tage of the fact that, for two electrons on an A-site lat- 
tice with on-site repulsion, it is sufficient to study the 
two-particle Green's function only in the A-dimensional 
Hilbcrt space of the doubly occupied sites, instead of 
the full A 2 -dimensional space 0. In the present pa- 
per we generalize such considerations for the dilute case, 
in which only / out of the N sites are replaced by impuri- 
ties. These impurities may have different single-electron 
(le) energies {e^} and on-site electron-electron (e-e) in- 
teractions {Ui}. As we show, the exact eigenenergies of 
the two electrons are found from the eigenvalues of a 
(small) I x I matrix, which involves the eigenfunctions 
and the eigenvalues of the le Hamiltonian. One surpris- 
ing result shows that in general the interacting eigenen- 
ergies are bounded between consecutive non-interacting 
two-electron (2e) energies, so that the energy cost due 
to the interactions is usually much smaller than the av- 
erage repulsion (U), asigned to each pair in the Coulomb 
blockade approach (?]]. 

Our study is also relevant for quantum dots ||. A 
quantum dot coupled to electrodes has been modeled by 
one impurity (I = 1) on a one dimensional (ID), iV— site 
chain M. This represents a special case of the Anderson 
model PI , with a momentum-dependent hybridization 



between the impurity and the conductance band. Unlike 
the case of a momentum-independent hybridization, we 
find that the behavior of the two electrons on the "dot" 
has a rich phase diagram, as function of the "dot" site 
energy eo and the hybridization, i.e., the hopping energy 
to between the "dot" and the leads. These parameters 
can be tuned experimentally, by varying the voltage on 
a gate coupled capacitively to the dot and the barriers 
between the dot and the leads Jll| . We also find a der- 
ealization transition at large U, for sufficiently negative 
eo. Our method allows a similar analytic treatment of 
the 1=1 case in higher dimensions. 

More complex examples include a cluster of I impuri- 
ties coupled to ID leads, representing one large dot p2[ , 
two separate impurities, representing double quantum 
dots (or artificial "molecules") |13|, etc. The latter case 
may also shed light on the nature of the low-lying states 
in the impurity band of doped semiconductors, which 
involve "molecules" of impurity-pairs [jl4|. For 1 = 2, 
the "molecule" adjusts to our bound on the interaction 
energy by crossing over to the Mott antiferromagnetic 
(singlet) state, in which the electrons are localized on 
separate ions |L4|,|l5). Finally, the case I = N repre- 
sents the 2e Hubbard model, whose non-random version 
(e, = 0, Ui = U) can be solved analytically in general 
dimensions using our method. 

We start with the le Hamiltonian 



Ho — 



E 

(n,m) 



t \n)(m\ + h.c.) + ^2ei 



(1) 



where \i) is a (spin-independent) state fully localized on 
site i, the first sum runs over all the site pairs in the sys- 
tem (including the impurities) and the second sum runs 
over the I impurities. We first find the le eigenstates 
| a) = ^2 n 4> a {n)\n) and eigenenergies e a of Tio- This is 
usually easy, being a linear problem (in any case, this 
involves at most the diagonalization of an N x N ma- 
trix). The on-site interactions take place only on the 
impurities, 



Hint = £f =1 f7in 



(2) 



where hi a = \i<r)(ia\ is the number operator of electrons 
with spin a =|, J, on site i. It is convenient to construct a 
basis for the 2e Hilbert space from the eigenstates of Ho- 
These are split into singlet (S) and triplet (T) spatial 
states, both with energies e a b = £a + £&: 
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|a6)s = (|a(l))|6(2)} + |a(2))|6(l)))C„ b! 

|a6) T = (|a(l))|6(2))-|o(2))|6(l)))/V2, (3) 

where \a(j)) represents the le eigenstate of Ho for elec- 
tron j, ( a b — 2~( 1+Sab ^ 2 (with the Kronecker <5), and 
the T states are used only for a ^ b. To avoid double 
counting, we consider only states with e a < e . 

Using this basis, the matrix elements of Hint involving 
T states vanish for the on-site interaction of Eq. (||) , and 
the energies of the T states remain e ao . We hence focus on 
the S states, and omit the subscript S. The nonzero S-S 
matrix elements are (ab\Hi n t\cd) = Sf =1 t r j77 a f,(i)*ry C( j(i), 
where r] a b(i) = 2£ oo </> o (i)0b(i). We now write the eigen- 
functions of Ho + H m t, with energy E, as Y, ea < eb x a b\ab). 
The coefficients must obey 

(E - e ab )x ab = Eiy/UiT)ab{i)* A, (4) 
where Ai = Z ec < ed VU~iVcd(i)x c d- Using @ we find 

Ai = T<j =1 SijAj, (5) 

where 

Sij(E) = yJWi ilab{i)Vab{])*l{E - c«t). (6) 

Sa/jTWi = G(jj;ii;E) = (jj\(E - Ho^U) is the 
non-interacting 2e Green's function in which the two 
electrons are on the same site, cf. ||. 

In addition to degenerate solutions |l6| , the new eigen- 
states are found by requiring a non-zero solution to Eq. 
(||) , namely that the / x / determinant D(E) = \ \ Sij — Sij \ \ 
vanishes. One way to find the new eigenvalues {E} is to 
find the / eigenvalues Si(E) of the matrix S, and then 
solve the equations Si(E) = 1. The A^s are then given 
by the eigenvectors related to Si(E), and the {x a b(E)Ys 
are found (up to a normalization constant) from Eq. (|J) . 
We emphasize that, unlike a perturbative expansion, this 
formalism gives exact values for E even for very large 
Ui's. As we show below, this enables us to find interest- 
ing transitions in that limit. 

The equation D(E) = can be used to find bounds 
on the energies {E}. Since each Sij has poles at ev- 
ery £a6, we generally expect each Si also to have such 
poles. Although the residues of some of these poles 
may vanish in special circumstances (e. g., Si = Su = 
Ui/(E — 2ei) for isolated impurities, t nm = 0, or for spe- 
cial symmetric cases, see below), this is not expected to 
happen in the general random case, which we discuss 
now. When E is very close to a non-interacting en- 
ergy e ab , one has Sij ~ y/U^r) ab (i)r] ab (j)* / (E - e ab ). 
It is then easy to show that in this approximation, 
D{E) » (-l)'(l - TrS) » (-1)'[1 - (U) ab /(E - e afc )], 
with (U)ab = J2i U i\Vab(i)\ 2 - Indeed, for very small £/, 
the equation D(E) — reproduces the lowest order per- 
turbation result, E w e ao + (U) a b- As E crosses through 



e a b, (~1) I+1 D(E) jumps from — oo to oo. As E in- 
creases between two consecutive non-interacting eigen- 
values, {-1) I+1 D(E) has no singularities, and thus varies 
smoothly from oo to — oo. Except for the special cases 
mentioned above, for which D(E) splits into products 
involving subspaces of the levels (see below), we thus 
conclude that D(E) = must have at least one (and 
up to /) solution(s) between every pair of such consecu- 
tive energies, and the new energies maintain the sequence 
of the non-interacting ones. Specifically, if the Ic low- 
est levels are e g and e u then the non-interacting lowest 
energies are e gg = 2e g and e gu — e g + e u , and the inter- 
acting ground state energy E gg obeys e gg < E gg < e gu . 
Thus, the effective interaction cost for adding the two 
electrons, A gg = E gg — e gg , is now bounded by {e u — e g ), 
and is always smaller than (U). For small (e u — e g ), this 
cost is negligible. Similar results apply to all the levels 
except the largest one (see below). This renormalization 
of the interaction energy has direct consequences for the 
issue of the Coulomb blockade, where one assumes that 
A a b ~ {U) a b, usually assumed to be independent of ab. It 
is interesting to note that small (and not evenly spaced) 
values of A were observed in a series of 2D quantum dots 
p7[ . It is tempting to relate these observations to our 
result. 

The simplest example concerns one impurity (or "dot" ) 
on a ID wire, closed with periodic boundary conditions. 
We place the impurity at site i — 0, with energy eo and 
with real matrix elements to,i — to,N—i = to- The other 
nearest neighbor matrix elements are set at t n ^ n +\ = t = 
1, for n = 1, 2, N — 2, and all energies are scaled by t. 
The eigenenergies of Ho have the form = 2cos(fc). 
Out of these, N/2 belong to odd states, of the form 
(j) k (n) oc sin(fcn), with k = 2irt/N, I = 1,2,..., N/2 (for 
convenience, N is even). These states do not "feel" the 
impurity (</>fc(0) = 0), and remain unchanged for all eo 
and U . The remaining N/2 states are even functions of 
n, with 4> k (n) = fc (O)t o cos[fc(A-2n)/2]/cos(fcA/2) for 
n 7^ 0, where |0fc(O)| is determined by J2 n l < / > fe( n )| 2 = 1- 
The allowed values of k are given by 

27tan(fcA/2) sin(fc) = 2(1 - 7) cos(fc) - e , (7) 

where 7 = (to A) 2 = ^ n ana lysis of this equation 
yields the surprising phase diagram shown by the full 
lines in Fig. 1: In region A all the states are delocalized 
(k is real). In region C (or D+F) there exists one bound 
state above (or below) the conduction band, with a lo- 
calization length I/k+ and energy e+ = 2cos1i(k + ) > 2, 
(or 1/k- and e_ = — 2cosIi(k_) < —2), where e K± = 
±e /2 + ^(e /2) 2 - I + 27. Finally, both bound states 
exist in region B. It is interesting to note that when eo 
is inside the original conduction band, |eo| < 2, the state 
on the impurity becomes delocalized for any infinitesimal 
hybridization 7, and the impurity does not imply local- 
ization (as might be anticipated in ID). Furthermore, 
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this state becomes localized for large 7. This surprising 
result arises since a larger 7 implies a larger repulsion of 
the localized state below the band. 

Combining the normalization condition with Eq. (0) 
yields (for TV > 1) |0 fo (O)| 2 = 1/[1 + 2 1 /{e 2K ~ 1)] for 
the bound states b = ±, and |<M0)| 2 = 2nf(k)/N, with 
7r/(fc) = 4 7 sin 2 (fc)/[47 2 sin 2 (fc) + (2(1 - 7) cos(fc) - e ) 2 ]. 
The latter result implies that unless one has a uniform 
chain, where 7=1 and 60 = 0, the weight of the states 
with k at the band edges and tt vanishes on the impu- 
rity. This will turn out to be crucial below. 

For I = 1, Eq. (^|) reduces to a single equation, 
D(E) = Us(E) — 1, and the eigenenergies E obey 
s(E) = 1/U. Similar equations are encountered in nu- 
merous cases, e. g., in the Cooper 2e problem or in 
Kohn's model of the insulating state [Ell. Indeed, s(E) 
jumps from —00 to 00 as E crosses each non-interacting 
energy e a b, and we find one new eigenvalue between every 
pair of non-interacting energies, as described above. In 
regions B, D and F, e g = e_ and e u = —2 (the bottom of 
the band). Thus, A gg is bound by |e_| — 2, and is always 



smaller than U. In fact, A„ 



when the le bound 



state approaches the band, at the line eo = 2( 7 — 1). For 
another insight into the smallness of A gg , rewrite s = 
\<f>-(0)\ 4 /(E-2e-)-T(E). For E < e gu = e_-2, one has 
T > 0, and |0_(O)| 4 /(.E-2e_) = l/U eS , with the smaller 
renormalized repulsive energy U e s = U/(1+TU) < U. It 
is interesting to note that in regions D and F, T is pro- 
portional to 7, and therefore the renormalization of U 
increases with the hybridization. The Coulomb blockade 
picture is restored for 7 <C 1. 

The details of s(E) depend on the parameters eo and 
7. To treat the bands (of states \kk') and | ± k)) in the 
limit N — > 00, we replace the summations by integrals: 



E 



2 1 -W|&(0)<M0)|< 



b.b 



E-e b -e' b 
2^|<M0)| 2 f dk 



f(k) 



E — 2 cos(fc) — e& 

"2 I dk f dk' — ' 

Jo J E-2 cos(fc) - 2 cos(fc') 



(8) 



where the bound states b, b' = ± are included only when 
they exist (regions C, D, F and B). For E inside the 
band we must return to the discrete sum, and the dense 
energies practically don't shift. 

Our most interesting results arise in region D+F. Here 
we have only one bound state, b = — , so that s(E) is 
negative for E < e gg = 2e_, and decreases from +00 
as E increases above this value, towards the band which 
starts at e gu = e_ — 2. For finite N, s would diverge 
towards —00 as e gu is approached, implying a persistent 
bound solution with a discrete energy below the band 
for all U. This would also be the case for the contin- 
uum case, N — > 00, if one had a non-zero value of /(0) 



(from |<^>fc(0)| 2 at k = 0), due to the divergence of the ID 
density of states there. However, as noted above, f(k) 
vanishes at k = 0, yielding a finite value s(e gu ) = s c . 
We studied s c as function of eo and 7, and found that in 
region D of Fig. 1 one has s c < 0, so that the equation 
s(E) = 1/U still has a discrete bound ("insulating") state 
there. However, in region F one has s c > 0, implying a 
disappearance of this bound state for U > U c = l/s c . 
This transition is intuitively easy to understand: when 7 
is very small, this transition occurs when the energy of 
the two electrons in the isolated atom, 2eo + U, exceeds 
that of the state in which one electron remains bounded 
while the other moves to the band, which is equal to 
eo — 2. At finite 7, a larger negative eo implies a smaller 
localization length The electrons are then more 

localized on the impurity, and A gg is larger, leading to 
the "insulator to metal" transition from region D to re- 
gion F. It should be noted that although the bound state 
is a singlet, with total spin zero, the new ground state 
in region F has one bound electron and one "free" elec- 
tron. Such a state does not feel the e-e repulsion, and is 
thus practically degenerate with the slightly lower triplet 
state (for large TV, the difference is of order l/N). Unlike 
the "insulator" singlet (or "antiferromagnetic" ) ground 
state, which has no net magnetic moment, this "metal- 
lic" state in region F is paramagnetic. This difference 
should be measurable in an external magnetic field. 

It is interesting to study the crossover from the "meso- 
scopic" case, of finite N, to the thermodynamic limit 
discusssed above. Firstly, for finite large N the tran- 
sition from region A to region D is smeared, occuring 
when (e_ + 2) becomes comparable to the spacing be- 
tween the band states, of order l/N 2 . Secondly, in the 
mesoscopic case, s(E) diverges to —00 as E approaches 
the lowest band state (slightly above e gu ). However, in 
region F, s(E) first gets very close to s c , and only then 
drops sharply to —00. For U > U c , the resulting "bound" 
state will thus occur very close to e gu , making it almost 
indistinguishable from the states inside the band. This 
implies a "smeared" transition from region D to region 
F for large mesoscopic systems. 

Similar interesting effects occur above the band, in re- 
gions A, D and F. The highest band energy is E = 4, 
corresponding to the upper limit of the last term in Eq. 
(||). Since /(tt) = 0, s(E) remains finite as E approaches 
4 from above, and < s(E = 4) = s x < 00 [s x — 00 only 
on the boundary eo = 2(1 —7)]. As E increases from 4 to 
00, s decreases from s x to zero, implying a discrete solu- 
tion to the equation s(E) = 1/U, with E > 4 above the 
band, if U > U x = l/s x . We found (by calculating x a b) 
that the corresponding wave function decays away from 
the impurity. Since this new state is above the band, it 
describes an "antibound" 2e excitation, which is remi- 
niscent of the upper Hubbard band at low concentration 
p9[ . Our analysis finds a phase transition of this excited 
state, from an insulating localized state at U > U x to 
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a conducting state for smaller U. This transition also 
becomes "smeared" for finite N. 

Now we give a very brief summary of our results for / = 
2, on the ID ring. Consider the very symmetric case, with 
the two impurities at sites i = i\, 12 having e; = eo, Ui = 
U , ti t i±i = to and 12 — i\ = R. The le real wave functions 
separate into even and odd ones, with 4> a (ii) = ±<X«2)- 
Depending on eo, to and R one again finds a rich le phase 
diagram. In particular, the low energy le states may 
start with the band, or have a bound even ("bonding") 
state or (above a minimal value of R, R x ) have also an 
odd bound ( "antibonding" ) state between the bonding 
state and the band. The high symmetry implies that 
5*11 = S22 and S12 — S21, with eigenvalues S±(E) = 
Sn ± Si2- Writing S±(E) as sums over non-interacting 
singlet states, we find that 5+ contains only the even- 
even and the odd-odd states, while <S_ contains only the 
even-odd states. Thus, D(E) = (S + — 1)(<S_ — 1), and the 
problem splits into two separate spaces. Consequently, 
A gg is now bounded by the first non-interacting excited 
state within the even-even and odd-odd subspace. The 
analysis of each equation S± — 1 is now similar to that 
of J = 1 , and the resulting phase diagram again contains 
derealization and magnetic transitions [^0| . It is possible 
to follow the coefficients {x a h.}, and see exactly how U 
causes a gradual crossover from the "molecular" to the 
"atomic" limits as R increases |l5| . In the latter localized 
states, the effective interaction energy obeys our bound 
by being of order 1/U. This crossover would be missed 
for R < R x , if one were to ignore the hybridization with 
the band (as commonly assumed pl[). 

Our method is easily generalized to higher dimensions 
and to larger /. It should be particularly helpful in nu- 
merical work on random systems. Another simple ex- 
ample concerns the non-random Hubbard model, where 
I = N. The matrix Sij is easily diagonalized in Fourier 
space, yielding the upper Hubbard band. This and other 
examples will be published elsewhere |2(J ■ 
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FIG. 1. eo — 7 phase diagram for the single impurity case. 
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